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Abstract 

A linked cluster expansion for the calculation of ground state observables of nuclei with realistic 
interactions has been developed. Using the V%' potential the ground state energy, density and 
momentum distribution of complex nuclei have been calculated and found to be in good agreement 
with the results obtained within the Fermi Hyper Netted Chain, and Variational Monte Carlo 
approaches. Using the same cluster expansion, with wave function and correlations parameters 
fixed from the calculation of the ground-state observables, various high energy scattering processes 
off complex nuclei have been calculated taking final state interaction effects into account by means 
of the Glauber multiple scattering series. 
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I. INTRODUCTION 

The knowledge of the nuclear wave function, in particular its most interesting and poorly 
known part, viz the correlated one, is not only a prerequisite for understanding the details 
of bound hadronic systems, but is becoming at present a necessary ingredient for a correct 
description of medium and high energy scattering processes off nuclear targets; these in fact 
represent a current way of investigating short range effects in nuclei as well as those QCD 
effects (e.^'.color transparency, hadronization, dense hadronic matter, etc) which manifest 
themselves in the nuclear medium. The problem of using a realistic many-body description 
is not trivial, for one has first to solve the many body problem and then has to find a way to 
apply it to scattering processes. The difficulty mainly arises because even if a reliable and 
manageable many-body description of the ground state is developed, the problem remains 
of the calculation of the final state. In the case of complex nuclei, much remains to be 
done in order to achieve a consistent treatment of Initial State Correlations (ISC) and final 
state interaction (FSI), which is feasible in the case of few-body systems while to date is still 
lacking for nuclei with mass number larger than A — A; the structure of the best trial wave 
function resulting from very sophisticated calculations {e.g.the variational Monte Carlo ones) 
is so complicated, that its use in the calculation of various processes at intermediate and 
high energies appears to be not easy task. We have developed an economic method for the 
calculation of the ground-state properties (energies, densities and momentum distributions) 
of complex nuclei within a framework which can be easily applied to the treatment of various 
scattering processes, keeping the basic features of ISC as predicted by the structure of 
reahstic Nucleon-Nucleon (NN) interactions. 

II. THE CLUSTER EXPANSION 

We write the nuclear Hamiltonian in the usual form, i.e.: 



where the vector x denotes the set of nucleonic degrees of freedom, and the two-body 
potential J2i<j V2{Xi,Xj) = Yln=i with r^j = |rj — rj\ denoting the relative 

distance of nucleons i and j, contains the well known spin- and isospin-dependent operator 




(1) 



2 



0^^\ The evaluation of the expectation value of the nuclear Hamiltonian is object of 
intensive activity which, in the last few years, has produced considerable results; our goal is 
to present an economical, but effective method for the calculation of the expectation value 
of any quantum mechanical operator A 

where 0o is a Shell-Model (SM), mean-field wave function, and F is a symmetrized correlation 
operator, which generates correlations into the mean field wave function. According to the 
two-body interaction in Eq. (Q), the correlated ground state wave function 0o is cast in the 
following form: 

A 

F{xi, X2 ... xa) 0o(a3i, X2 ...xa) = S Yl fiXij) 0o(a:i, X2 ... xa) (3) 

with 

TABLE I: Benchmark calculation for the ground state potential, kinetic, total energy and energy 
per particle for ^^O, at the first order of r/-expansion (Eq. ©), compared with the FHNC results 
for the corresponding quantities calculated with the same wave function. 
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-390.37 323.50 


-66.87 -4.18 




FHNC f5] 


-390.30 325.18 


-65.12 -4.07 
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In the following, we are j 
we used to evaluate Eq. ( 


^oing to describ 
21); the solution 


e the cluster expansion (77 — expansion) technique 
can be found by applying the variational principle. 



with the variational parameters contained both in the correlation functions and in the mean 
field single particle wave functions. The expectation value (A) defined in (121) can be cluster- 

nnn 

expanded; at first order of the 77— expansion I2, it reads as follows 

= (0o I E {k^k - ^) I 0o) - (i)o(<^o I E {f^lk - 1) I <Po) , (5) 
i<j i<j 

where {A)o is given by {(j)o \A\ (po)- The 2nd order term can straightforwardly be obtained 
by the same technique used to derive Eq. (0). Given the two-body interaction as in Eq. ((TJ, 



the expectation value of the Hamiltonian can be written as: 

Eo= f dk^nik) + E / dr,dr2V^^\n2)p^^){rur2), (6) 

where T2) is the (spin and isospin dependent-) Two Body Density (TBD) matrix [l|, 

Ql, and n{k) is the nucleon momentum distribution, defined in terms of the non-diagonal 
One Body Density Matrix (OBDM), p^'^\ri,r[), as 

n{k) = ^ / dridr;e^^-(^-^'^) p{r^,r',) . (7) 

The knowledge of the OBDM allows one to calculate other relevant quantities like e.g. the 
density distribution p{r) = p^^\ri = r[ = r). The results of calculation of the ground state 
energy, the charge density and the two-body density and momentum distribution using the 
reahstic V8' interaction P| is discussed in detail in Refs. [ll, O, S] ; in Table HI some results 
for a benchmark calculation of the quantity in Eq. ^ are shown and it can be seen that the 
agreement with the reference FHNC calculation is very good. Results for the momentum 
distribution and radial TBD, defined as /5(„)(r) = / dR p\^^'^{r_\,r2), where r and R are 
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the relative and center of mass vectors of particles 1 and 2 [jj, |2[, are shown in Fig. 
we obtained the best wave functions for ^^O and ^°Ca, used to calculate the quantities 
shown in Figs. ^ and |21 and to be used in various applications, from the results achieved 
variationally within the FHNC jsl approach, i.e., we used the correlation functions provided 
by the FHNC method and then optimized the single particle wave fucntions in order to 
reproduce the experimentally observed charge densities and radii for these two nuclei; we 
found that ^^O can be described both with with harmonic oscillator or Saxon- Woods mean 
field functions while, in the case of ^°Ca, Saxon- Woods functions have to be used to produce 
satisfactory densities and momentum distributions, though the agreement with the FHNC 
results for the ground state energy is not as satisfactory as in the ^^O case . The method 
we have developed appears to be a very effective, transparent and parameter-free one. An 
extensive discussion of the applications of the method of the cluster expansion to scattering 
problems is out of the scope of the present report and will be presented elsewhere. In this 
report, we show what the effect of the Glauber final state interaction on the distorted 
nucleon momentum distribution which appear in the semi-inclusive A{e, e'p)X reaction; i.e. 

"^(P-) = / e^^'"-^^^-^'^VD(ri,r;)rfricir;, (8) 




FIG. 1: The momentum distribution (Eq. ((71); left) and the radial two-body density (Eq. 
right) of ^^O calculated within the cluster expansion. The ry-expansion (thick solid) result for the 
momentum ditribution is compared with the VMC (squares) and the FHNC _5] (stars) ones; 
the mean field (dots) and the Jastrow (thin solid) results is also show, in the left figure. The six 
components of p^^-^ 3]) which couple with the first six components of the ^8' potential are shown 
in the right figure. 



where = q — p is the missing momentum, defined as the difference between the three- 
momentum transfer and the momentum of the detected proton, and 

{^o\S^p{ri,r[)S'\^oy 



(9) 
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FIG. 2: The distorted momentum distribution of Eq. Q for ^^O (left) and ^"Ca (right). Thick 
solid line: no final state interaction taken into account; dotted, dashed and thin solid lines: FSI 
in the Glauber framework for perpendicular, parallel and antiparallel kinematics, respectively; 9 
represents the angle between the missing momentum p„ and the three-momentum transfer, q. 



5 



is the distorted OBDM. The FSI is introduced by the Glauber approximation, i.e. where 
pi{fi,f[) is the OBDM operator, S is the Glauber S-matrix, i.e. 

A A 

S ^ Sciri . . .r^) = n G'(ri,r,) = Y[[l - 9{z, - z^)r{h, - b,)] (10) 

j=2 3=2 

with asymptotic values of the parameters appearing in the profile F functions The 
distorted momentum distributions (jH} are presented in Fig. |21 

III. SUMMARY 

To sum up, we have shown that, using realistic models of the NN interaction, a proper 
approach based on cluster expansion techniques can produce reliable approximations for 
those diagonal and non diagonal density matrices which appear in various medium and high 
energy scattering processes off nuclei, so that the role of nuclear effects in these processes can 
be reliably estimated without using free parameters to be fitted to the data. Our approach 
has been applied to semi-inclusive A{e, e'p)X reaction ; we also calculated the total 
neutron-nucleus cross section P, Q| within the same Glauber approach, obtaining, once 
inelastic shadowing effects are taken into account, a very good agreement with data over a 
wide range of mass number A, which confirms NN correlations can play an inmportant role 
in high-energy scattering processes. 
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